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We propose that an excitonic gap can be generated along nodal directions by Coulomb interaction 
in the mixed state of d-wave cuprate superconductors. In a superconductor, the Coulomb interaction 
usually can not generate any fermion gap since its strength is weakened by superfluidity. It becomes 
stronger as superfluid density is suppressed by external magnetic field, and is able to generate a 
gap for initially gapless nodal quasiparticles beyond some critical fleld He- By solving the gap 
equation, it is found that the nodal gap increases with growing field H, which leads to a suppression 
of thermal conductivity at zero temperature. This mechanism naturally produces the field-induced 
thermal metal-insulator transition observed in transport experiments. 
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I. INTRODUCTION 

The low-energy spectral and transport properties of 
nodal quasiparticles in d-wave cuprate superconductor 
are very important issues. In the absence of external 
magnetic field, the ground state is occupied by uniform 
superconductivity. For a clean superconductor, the low- 
energy density of states vanishes linearly as N{uo) ~ |a;| 
upon approaching the Fermi surface. It acquires a finite 
value at zero-energy in the presence of weak impurity 
scattering. In this case, the nodal quasiparticles exhibit 
universal transport behavior in the sense that the electric, 
thermal and spin conductivities are independent of im- 
purity scattering rate at the w ^ 0, T ^ limii^ii^i^. The 
zero temperature thermal conductivity k is of particular 
interests since it is not affected by vertex corrections^.. 
Remarkably, the predicted universal thermal conductiv- 
ity has been confirmed by heat transport measurements 
at optimal doping^i^. 

When placed in an external perpendicular magnetic 
field, the superconductor enters into the mixed state in 
the range Hd < H < Hc2- Inside the vortex cores, the 
superfluid current is significantly reduced by the mag- 
netic field. The low-energy fermionic excitations in the 
mixed state are expected to have rather different low- 
energy behaviors comparing with those in the uniform 
zero-field condensate^. As revealed by heat transport 
measurements, the thermal conductivity loses its univer- 
sality and depends on the impurity scattering rate^. In 
addition, on the underdoping side, it decreases as the 
magnetic field grows wg^^. There seems to be a field- 
induced thermal mctal-to-insulator transition in some 
underdoped cuprate supcrconductors^ii. These experi- 
mental results can be intuitively understood by assum- 
ing that the nodal fermions acquire finite mass gap in 
the mixed state. A phenomenological expression for the 
nodal gap was proposed^ to understand the field-induced 
reduction of thermal conductivity. However, the dynamic 
origin for the gap generation has not been discussed. 

The goal of this paper is to suggest a mechanism for 
opening the field-induced gap for the initially gapless 
nodal quasiparticles. Generally, this mechanism would 



be realized by some kind of fermion self-interaction or bo- 
son mediated interaction. Such interaction should have 
the following two features: it is weak enough to be ir- 
relevant in the uniform superconducting state; it gets 
stronger with growing perpendicular magnetic field so 
that a finite gap is generated beyond some critical mag- 
netic field He. 

Qualitatively, the U(l) gauge fluctuation arising from 
strong correlation provides a good candidate for such 
mechanism. The t-J model of cuprate superconductors 
can be theoretically treated by the slave boson method. 
After making mean-field analysis and including fluctua- 
tions, there appears an emergent U(l) gauge field which 
interacts strongly with spin-carrying spinous and charge- 
carrying holons^. The superconductivity is realized by 
holon condensation below Tc, while the d-wave energy 
gap is formed by spinon pairing. In the superconducting 
state, the low-energy elementary excitations are gapless 
nodal spinous and the U(l) gauge boson is gapped via 
the Anderson-Higgs mechanism. The finite gauge boson 
gap weakens gauge interaction, so usually no fermion gap 
can be generated. However, once the superfluid density 
is suppressed by external magnetic field, the gap of U(l) 
gauge boson decreases and the strength of gauge inter- 
action increases with growing magnetic field^. Then a 
finite gap for nodal spinons could be generated by gauge 
fiuctuation, leading to suppression of thermal conductiv- 
ity. Unfortunately, it is hard to average over the vortex 
distributions within this formalism due to the complexity 
brought by spin-charge separation. 

The gapless nodal fermions might acquire a gap via the 
magnetic catalyst mechanisntii when they are placed in 
an external magnetic field. But this mechanism depends 
on a crucial assumption that the fermion stays in the 
lowest Landau leve l^^'^^ . However, in the case of high 
temperature superconductor, the Landau level has been 
shown not to be the appropriate description of fermion 
energy spectrum in the mixed statei^. Therefore, the 
magnetic catalyst mechanism is unlikely to be at work. 

In this paper, we study the possibility of gap gen- 
eration due to the long-range Coulomb interaction be- 
tween charged nodal quasiparticles. Two quasiparticles 
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that carry the same charges always experience a repul- 
sive Coulomb force, while the quasiparticle and quasihole 
experience an attractive Coulomb force. When the at- 
tractive force is sufficiently strong, it is possible that a 
Dirac quasiparticle is combined with a Dirac quasihole to 
form a stable excitonic pair. Through this mechanism, 
the gaplcss fermion acquires a finite excitonic gap. 

Recently, this kind of gap generation was argued 
to lead to an insulating ground state in single layer 
graphene, when the Coulomb interaction strength g is 
larger than some threshold gc and the fcrmion flavor N 
is less some threshold A^ ^^^i^^i^^ . Moreover, an interesting 
superfluidity was predicted to exist in bi-layer graphene 
based on a similar paring instability—. The long-range 
Coulomb interaction is also very important in cuprate 
superconductors. First of all, it lifts the gapless Gold- 
stone mode up to plasmon mode, which is actually the 
rudiment of Higgs mechanism. Its importance in the for- 
mation of stripe phase has been emphasized by several 
authorsi^. However, its role and influence on nodal quasi- 
particles arc still in debate^. In the absence of a reliable 
microscopic theory of Coulomb interaction, we resort to 
the phcnomcnological approach. 

From the available extensive experiments, we know 
that the nodal quasiparticles have rather long mean free 
path and behave like well-deflned Bogoliubov-Landau 
quasiparticles in the uniform superconducting state^i^i. 
This fact and its excellent agreement with BCS-type 
analysisii^i^S implies that the Coulomb interaction must 
be fairly weak in the superconducting state and gener- 
ally can not generate any fermion gap, except in the 
lightly doping region. On the other hand, in the non- 
superconducting ground state, it is generally believed 
that there are no well-defined Landau quasiparticles. It 
is reasonable to expect that the long-range Coulomb in- 
teraction is very strong in this state. The field-induced 
mixed state lies between these two extreme limiting 
cases. As the superfluid density decreases with magnetic 
field, the effective strength of Coulomb interaction gets 
stronger. For sufficiently strong interaction, a dynami- 
cal fermion gap can be generated by forming excitonic 
pairs. To implement this intuitive picture with ( explicit 
) computations, we assume a phenomenological form for 
the effective interaction strength which is a function of 
magnetic field H. After solving the associated gap equa- 
tion, we find that the growing magnetic field drives the 
system towards a phase transition into excitonic insulat- 
ing state beyond some critical value He- Once the nodal 
fermion acquires a finite gap m, the low-energy fermionic 
excitations are significantly suppressed below the scale 
m, leading to reduction of thermal conductivity. ( This 
can help to understand the transport behaviors observed 
in the mixed state of cuprate superconductors. ) 

Besides the thermal metal-insulator transition, an- 
other important issue about field-induced phenomena 
is the enhancement of antifcrromagnctic correlations 
inside the vortex cores. Such microscopic coexis- 
tence of magnetic order and superconductivity has 



been investigated experimentall3 r^^i^'^i^^i^^i^^i^'i^^ and 
theoretically22i^. From the field theoretic point of 
view, the field-induced antiferromagnetism or spin den- 
sity wave can be represented by a mass term for nodal 
fermion a^'^^ . With this identification, the mechanism re- 
sponsible for the thermal metal-insulator transition can 
also account for the existence of field-induced magnetic 
order in the mixed state. 

Before presenting the technical details, we would like to 
point out that a number of assumptions and approxima- 
tions will be used to simplify discussions on the issue of 
dynamical gap generation. Thus, the conclusions reached 
in this paper are reliable only at the qualitative, rather 
than quantitative, level. 

The paper is arranged as follows. In Sec. II. we build 
the model and write down the gap equation. In Sec. 
Ill, we propose the phenomenological form of the effec- 
tive Coulomb interaction and calculate its dependence 
on magnetic field H. In Sec. IV, we solve the gap equa- 
tion and give the field dependence of critical coupling 
gc and dynamical gap. The qualitative understanding of 
transport experiments is also discussed. We ends with a 
summary and discussion in Sec. V. 

II. MODEL HAMILTONIAN AND GAP 
EQUATION 

We begin our discussion with the following Hamilto- 
nian of d-wave superconductor 

^0 = 5Z*k[^kT3-AkTi]$k, (1) 
k 

where the standard two-component Nambu spinor rep- 
resentation = (cj^,c_kx) is adopted and is PauH 
matrix. The electron dispersion is Ck = —2t[coskxa + 
coskyo) — with /.io being the chemical potential and 
the o?-wave energy gap is Ak = ^{cos kxCi — cos/cyo). 

The quasiparticle spectrum is i?k = i/e^ -I- A^, which 
has four nodal points at the Fermi level. Linearizing the 
dispersion in the vicinity of th e nodes, one obtains the 
spectrum i?k = a/ Vpkf -\- "u^fcl, where ki is perpendicu- 
lar to the Fermi surface and k2 is parallel to the Fermi 
surface. The four-component Dirac spinor can be defined 

*I(2)(q>w„) = (c|(k,w„),c^(-k,-w„), 

cj(k - Qi(2), w„), cx(-k + Qi(2), 

where Qi(2) = 2Ki(2) is the wave vector that connects 
the nodes within the diagonal pairs, k = Ki-|-q with q ^ 
Kj. Here, we use the four-component spinor because it is 
impossible to define chiral symmetry in two-component 
representation of fcrmion field in (2+1) dimensions. 

The continuum Hamiltonian of free Dirac fcrmions can 
be written as 

72WAi9i/)^i + (1 ^ 2, a; 2/)i(2) 
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where ^ = 4" ^70- The 4x4 matrices can be chosen as 
7o = CTi (g) To, 7i = —ia2 ® T3, and 72 = ia2 <8) n, where 
ai acts in the subspace of the nodes in a diagonal pair, 
Ti acts on indices inside a Nambu field. There are two 
matrices anticommuting with them, 73 = ia2 (8) T2, and 
75 = (73 <E) tq. The matrices satisfy the Dirac algebra 
{7^,7,} =2diag(l,-l,-l). 

The Hamiltonian for the Coulomb interaction is 



replace the vertex function Tq by 70. Now the complete 
propagator can be formally written as 



He 



1 ^ 
Att ^ 



*,(r)7ovE',(r)^— i',,(r')7o^'.'(r'), 



r — r' 



(3) 

The bare Coulomb interaction in momentum space is 
^(9) = The parameter g measures the strength 

of bare Coulomb interaction in the non-superconducting 
ground stated. It is easy to see that this model is 
very similar to that in single layer graphene, where the 
low-energy excitations are also massless Dirac fermions. 
Unlike the semi-metal background in graphene, in the 
present case the Dirac fermions are in a charge conden- 
sate, which unavoidably affect their properties and the 
Coulomb interaction between them. 

From the results obtained in graphene, it is known that 
sufficiently strong Coulomb interaction can lead to finite 
excitonic gap, when the fcrmion flavor is below a critical 
valuei^ii^ii^. We speculate that a similar pairing instabil- 
ity to occur in the mixed state of d-wave superconductor. 
When studying the gap generation, we fix the physical 
fermion flavor N = 2. Thus, coupling g becomes the 
only variable that tunes the excitonic phase transition. 

The Hamiltonian is invariant under the continuous chi- 



ral transformation 5* 



^'^i. It will be dynamically 



broken once the Dirac fermion acquires a finite mass via 
the effective Coulomb interaction. This phenomenon is 
non-perturbative in nature and generally can be studied 
by analyzing the self-consistent Dyson equation 

G-\p)^G^\p)+T J ^,Gik)To{p,k)Vip-k), 
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(4) 

where the (2-f l)-dimensional momentum is defined as 

fc=:(iLi;„,k). (5) 

The Matsubara frequency is tOn = {2n+l)TTT for fermions 
and LUn = 2n7rT for bosons. Here, ro(p, fc) is the full 
vertex function. The free propagator for massless Dirac 
fcrmion is 



Goik) 



1 



iiUnJO - Vpkiji - VAk2j2 



(6) 



Due to the Coulomb interaction, it becomes the com- 
plete propagator G{p), which is determined by the Dyson 
equation. In the case of QED3, the dynamical chiral sym- 
metry breaking can be most conveniently studied using 
the expansion^i. Here we follow the same strategy 
and keep only the leading order of 1/N expansion. So 
we can neglect the wave function renormalization and 



G(fc) 



1 



iuJnlo - vpkiji - i'Afc272 - m{k) ' 



(7) 



where m{k) denotes the Dirac fermion mass. Further, as 
shown in the context of QED3, at least at low energies 
and to the leading order of 1/A^ expansion, the velocity 
anisotropy is irrelevant to the critical behavior^. We 
simply set vp = va = ^ whenever they multiply the 
momenta (fci, ^2) in the gap equation. 
The full Coulomb interaction function is 
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ci\/g + Nx{q)' 



(8) 



The polarization function xil) contains all informa- 
tion about how the Dirac fermions response to the 
many-particle system. We first consider the non- 
superconducting ground state. Within the random phase 
approximation, the fermion propagator G and the vertex 
fimction Tq are both replaced by the bare ones, i.e.. 



Tr[7oGo(fc)7oGo(fc-(7)]. (9) 



Inserting the expression for the interaction function, the 
gap equation can now be written as 



m(fc) 



1 



tY Lol + 'k^ + m?{k)\(\\lg + NxiqY 

(10) 

where q = p — k. In the instantaneous approximation, 
xil) has the following zero frequency expression^ 



2T 

x(q) = — / dx log [2 cosh 



^a;(l - a;)|q|- 



2T 



(11) 



Now the gap m is independent of frequency and the fre- 
quency summation can be carried out with the result 



reft „(r,.) ..nh^^g^ 

(12) 

In the limit of zero temperature, x(q) — Iq|/8, the gap 
equation further simplifies to 



m(p) 



d^k 



1 



B7r2 ^k2+,7i2(k) |q|/5 + ^|q|/8 ' 



(13) 



The nontrivial solution m{p) of this integral equation 
signals the occurrence of dynamical mass generation. 

When the ground state is occupied by the uniform su- 
perconductivity, the Coulomb interaction function must 
be modified. In the mixed state, the superfluid density 
is a non-uniform quantity which has different values at 
different spatial positions. To study the gap equation in 
the mixed state, we should average over the vortices and 
obtain a mean value of superfluid density. This is the 
task of the next section. 
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III. ANSATZ OF EFFECTIVE INTERACTION 

In the underdoping and optimal doping regions, the 
ground state is occupied by the superconductivity, which 
significantly weakens the Coulomb interaction between 
Dirac fermions. Such effect can be described by calculat- 
ing the polarization function that incorporates the effect 
of finite superfluid density. However, it is not clear how 
to correctly calculate the polarization function in the su- 
perconducting state, so we will assume a phenomenolog- 
ical form for the effective interaction function. From the 
experimental facts, we know that the interaction must 
reaches its minimal value when the superfluid density 
takes its maximal value. As Ag decreases with growing 
magnetic field H, the effective interaction strength in- 
creases and eventually takes its maximal value after the 
superconductivity is completely destroyed. We assume 
the following ansatz for the effective interaction function 
in the superconducting state 



V{q,H) 



(14) 



where a is an adjustable parameter. This is the simplest 
function that can describe the reduction of strength by 
superfluid density. It is largest at the limit A^ = 0, 
and is smallest at the limit H = 0. The function 
(1 -I- NVo{q)x{q)){i + ct-As{H)) can be considered as the 
effective dielectric function of superconducting state. In 
the mixed state, the field-dependent superfluid density 
As(iJ) controls the effective strength of Coulomb inter- 
action. The parameter a must be properly chosen so that 
a moderately strong magnetic field He separates the gap- 
less and gaped phases. In order to see how the critical 
point depends on H, we need to solve the gap equations 
after including As{H). 

To study the gap equation, the superfiuid density 
As{H) should be obtained by averaging over the vortices. 
In the mixed state, the low-energy properties of d-wave 
superconductor are dominated by the extended quasipar- 
ticles in the bulk material, unlike the case of conventional 
s-wave superconductor. Volovik^ proposed a scmiclassi- 
cal approach and showed that the density of states varies 
as ^/H at low temperatures, which has been observed by 
experiments^. Within the semiclassical treatment, the 
effects of circulating supercurrent around vortices can be 
represented by a Doppler shift^i^i^ in the quasipar- 
tilce spectrum, oj a; -I- k • Vs(r), where Vs(r) is the 
superfiuid velocity at a position r and k is the quasi- 
particle momentum which can be approximated by its 
value at the node. Then the fermion Green function 
can be written as G'((jj,k, r) = G{lli + ei(r),k). where 
ei(r) = ki • Vs(r). The local value F{r) of any phys- 
ical quantity F determined by the Green function can 
be obtained using the above local Green function. The 
field-dependent quantity F{H) is written as the following 
spatial average F{H) = -j J (PvF{uj + e(r)), where the 
integral is taken over a unit cell of the vortex lattice with 
area A. Such averaging integral depends on the vortex 



distribution. The field-dependent quantity is 

/oo 
deF{uj + e)Vie), (15) 
-oo 

with probability function '^(e) = ^ / (PrS{e — k • V5(r)). 

There are several possible choices of V{e), which were 
discussed in Ref.-^. For example, the distribution func- 

tion of vortex liquid or solid is V{e) ~ 2(e^+E'^ ' 
for disordered vortex state, it takes the form 7^(e) = 

_..v 



exp 



The typical energy scale of Doppler 



shift \&Eh = ^ = f'vC'' ^^"^^"^ ^ " ($o/7ri/)^/2 is 
the radius of the unit cell of vortex lattice and $o — hc/2e 
is one quantum of magnetic fiux. The field dependence 
of a physical quantity, such as density of state or spe- 
cific heat, depend somewhat on the choices of distribu- 
tion function, but the qualitative result is not sensitive 
to the choice. 

The computation of superfluid density As{H) within 
the semi-classical approximation has already been per- 
formed in Refi^, so we just list the basic steps and cite 
the results. The superfluid stiffness is given by^i^^ 



Al^iT,H) = T'^ -A]^{T,H), 



(16) 



where r,;j is the diamagnetic tensor and A„ represents the 
normal fluid density divided by the carrier mass. In the 
Matsubara formalism, the normal fluid density is^ 



oo „ 
„ — — , c'HBZ 



t'F(k)t'Mk) 



Tr[G'(iw„,k)7o75G(iw„,k)7o75], (17) 

with the i-component velocity v'^p. Here, HBZ means the 
halved Brillouin zone of the d-wave superconductors, i.e., 
the domain with two neighboring nodes. The averaged 
normal fluid density is given bj*^ 



Kim = 



deV{e) 



duj tanh ^ ^ 

2T 4i7r 



HBZ (2^)' 

xTr[GA(t^ - e,k)7o75GA(w - e,k)7o75 
-Gr{uj - e,k)7o75Gfl(w - e,k)7o75]. (18) 



In principle, the superfluid density should be calcu- 
lated by including the complete fermion propagators 
Gr^a{'^, k) into (18). Hence it actually satisfles an equa- 
tion that couples self-consistently to the gap equation. It 
is not an easy task to solve these coupled equations in 
practice. However, there is a remarkable simplification if 
we are mainly interested in what happens in the vicinity 
of the critical point of chiral phase transition. Near the 
bifurcation point, the gap equation can be linearized and 
the gap appearing in the superfluid density can be taken 
to be zero. In the limit m ^ 0, the normal fluid density 
finally becomes^^ 



An{H) = 



VF 



deV{e)J{m,e), 



(19) 
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with J(e) = 2|e|. If we adopt the distribution function 
T'{e) of vortex Hquid, then the superfluid density 



K{H)=T- 



TTVA 



E 



(20) 



Here, r is the zero temperature superfluid density in the 
absence of magnetic field and the ratio between vp and 
f A appears as a coefficient. Although the anisotropy is 
irrelevant when uj;^(a) multiplies particle momenta, the 
ratio might be important in this expression. Wc simply 
set the ratio vf/va = 20 in the following discussions. 



IV. FIELD DEPENDENCE OF CRITICAL 
COUPLING AND THERMAL CONDUCTIVITY 

After getting the effective Coulomb interaction, we 
now study the self-consistent gap equation 
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FIG. 1: The dependence of critical strength gc on magnetic 
field H for several choices of a. 



77l(p) = 



d^k 



7Tl(k) 



1 



I 



87r2 + m2(k) |q|/5 + iV|q|/8 1 + aA,(i?) 



This integral equation can be solved numerically by the 
parameter embedding method, with the Coulomb inter- 
action strength g being the turning parameter. There 
is a critical value gc that separates the chiral symmetric 
phase (to — for g < gc) from the symmetry breaking 
phase (m ^ for g > gc). gc is just the critical point of 
chiral phase transition. We can see that excitonic pairing 
is quite different from conventional BCS type pairing: the 
former is produced only by sufficiently strong, attractive 
Coulomb force between particles and holes, while the lat- 
ter is triggered by arbitrary weak attractive force between 
electrons. The excitonic gap breaks the chiral symmetry 
and leads to the formation of antiferromagnetism or spin 
density wave in the field-induced vortex state. 

Following Refi^, the zero temperature superfluid den- 
sity is taken to be r = 1500K and the energy scale 
Eh ^ 30\/i?KT^^/2^ Xhe parameter a is a variable (in 
unit of eV^^) depending on doping concentration and the 
type of cuprate superconductors. It surely is not a uni- 
versal quantity and can not be uniquely determined. For 
completeness, we consider a number of possible values, 
a = 4, 5, 6, 7. For each value of a, the relationship be- 
tween critical strength gc and magnetic field H is shown 
in Fig. 1. As it turns out, critical strength gc decreases 
as the magnetic field H grows up. 

To judge whether an excitonic gap is generated for 
nodal fcrmions, we can simply compare the physical 
strength g with the critical value gc- Admittedly 
the exact value of physical strength g in the non- 
superconducting ground state is unknown. However, 
we can make a simple comparison between its value 
in cuprate superconductors with that in single layer 
graphenc. In graphcne. the typical value g is about 20. 
It is not unreasonable to estimate that the parameter g 
in the non-superconducting ground state be larger than 



20, since the correlation is known to be very strong in 
cuprate superconductors. 

If we assume that g = 50, then the critical field 
{21)Hc = 0,1.0,3.5,6.0 for parameters a = 4,5,6,7 respec- 
tively. For a = 4, a finite excitonic gap is generated even 
in zero field case. This is able to explain several experi- 
mental results performed in some underdoped cupratcs in 
the absence of external field H: the finite nodal gap found 
by photoemission^, the suppression of thermal conduc- 
tivity K from the universal value with lowering doping 
concentration^, and the coexistence of competing mag- 
netic order with superconductivity^. For larger values 
a = 5,6, 7, the gap is generated only for magnetic field 
H larger than some critical value He- These parame- 
ters are relevant to the doping regions in which the nodal 
gap and competing order appear only in the field-induced 
state22i2^i^i2^i2^. It appears that the phenomenological 
parameter a should be an increasing function of doping 
concentration. The cases for other values of g and a can 
be analyzed similarly. 

The generated gap will surely affect all observable 
physical quantities, such as specific heat, electric and 
thermal conductivity. Here we are primarily interested 
in the zero temperature thermal conductivity k. If we 
assume a constant fermion gap m and a small impurity 
scattering rate Fimp, then the zero temperature thermal 
conductivity has the expression^ 



K r2 

imp 
imp ' 



(22) 



It is no longer universal and depends explicitly on im- 
purity scattering rate as well as on gap to. Obviously, 
in the weak impurity limit. Fimp <C to, the thermal con- 
ductivity is rapidly suppressed from its universal value 
by TO. To see how k varies with magnetic field H, we 
need to know the field dependence of gap m{H). To this 
end, we solved the gap equation (21) numerically^ and 
presented the results for a = 4 and g ~ 200 in Fig. 2. It 
is evident that the gap m{H) is an increasing function of 
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magnetic field H . Thus in the symmetry broken phase, 
as magnetic field H grows, the thermal conductivity is 
suppressed and the system undergoes a phase transition 
from thermal metal to thermal insulator at high field 
limit, which is qualitatively in agreement with transport 
experiments^!^. 
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FIG. 2: The dependence of the generated gap on momentum 
and magnetic field for a = 4 and g = 200. A is the momentum 
cutoff. 

As revealed by transport experiments, the field induced 
reduction of thermal conductivity occurs only in under- 
doped and optimally doped cuprate superconductors^. 
In the overdoped region, the thermal conductivity is 
found to increase with growing mag netic field H^. This 
can be understood by assuming that the d-wave super- 
conductivity responses differently to external magnetic 
field in underdoped and overdoped cuprates. On the 
underdoping side, the magnetic field only reduces the 
superfluid density inside vortex cores but leaves the d- 
wave energy gap essentially unchanged. Due to the ad- 
ditional excitonic gap along nodal directions generated 
by Coulomb interaction, the low energy nodal quasipar- 
ticles arc significantly suppressed, thus reducing the ther- 
mal conductivity. However, on the overdoping side, the 
magnetic field destroys the c?-wave energy gap by directly 
breaking the Cooper pairs. As a result, extended quasi- 
particles are excited from the condensate and the thermal 
conductivity increases with growing field H . 

V. SUMMARY AND DISCUSSION 

In summary, we proposed a mechanism to explain the 
field-induced reduction of thermal conductivity in the 
mixed state of d-wave cuprate superconductor. In this 
mechanism, a finite gap for nodal fermions is generated 
by the strong Coulomb interaction between gapless par- 



ticle and hole excitations. Since the Coulomb interaction 
is usually weak in the superconducting state, such gap 
generation is possible only after the superfluid density is 
reduced by strong external magnetic field. The excitonic 
gap reduces the thermal conductivity at low temperature, 
so there is a thermal metal-insulator transition driven by 
magnetic field. 

There are several effects that might change the criti- 
cal behavior of chiral phase transition. For example, the 
long-range Coulomb interaction can be screened by the 
finite zero-energy density of states produced by disorder 
scattering and/or vortex scattering. Such screening ef- 
fect reduces the possibility of fermion gap generationi^. 
Only in clean superconductor and at magnetic field much 
lower than the up critical field iJc2, this effect can be ig- 
nored. On the other hand, the gap generation can also 
be promoted by other mechanisms. If there are strong 
contact interaction between nodal fermions, the possi- 
bility of pairing instability is significantly enhanced due 
to the positive contribution from contact interactioni^. 
These competing effects can be included into the above 
calculations along the steps presented in Ref J^. 

Throughout the present paper, the thermal fluctuation 
effect is totally omitted. This effect actually plays at least 
three important roles. First, the thermal fluctuation ef- 
fectively excites quasiparticles out of the condensate and 
hence reduces the superfluid density rapidly. Secondly, 
these thermally excited quasiparticles lead to screening 
of the Coulomb interaction. In addition, the excitonic 
paring will surely be suppressed by thermal fluctuations. 
These effects compete with each other, making the situ- 
ation rather complicated. This is why we simply neglect 
the thermal effects and consider only nearly zero temper- 
ature. To make an extension to finite temperatures, all 
these three effects should be carefully analyzed. 

Finally, we must admit that in the present work we 
utilized a number of assumptions and approximations 
when studying the dynamical gap generation for Dirac 
fermions. The results obtained in this paper are only 
qualitatively reliable. In particular, we have not arrived 
at a quantitative determination of the critical magnetic 
field He- Generally speaking, He must be a function of 
doping concentration, temperature, and type of super- 
conductor sample. Unfortunately, it is difficult to quanti- 
tatively include any of these effects. We wish the present 
work will be put on a firmer theoretical ground in the 
future. 
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